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Abstract 

The non-equilibrium Green's function formalism for infinitely extended reservoirs coupled to a 
finite system can be derived by solving the equations of motion for a tight-binding Hamiltonian. 
While this approach gives the correct density for the continuum states, we find that it does not lead, 
in the absence of any additional mechanisms for equilibration, to a unique expression for the density 
matrix of any bound states which may be present. Introducing some auxiliary reservoirs which are 
very weakly coupled to the system leads to a density matrix which is unique in the equilibrium 
situation, but which depends on the details of the auxiliary reservoirs in the non-equilibrium case. 
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I. INTRODUCTION 



Electronic transport in mesoscopic systems has been studied intensively for several 
years^"^. For one-dimensional systems, the Landauer-Biittiker formalism has played a piv- 
otal role in this subject^. For a wire in which only one channel is available to the electrons 
and the transport is ballistic (i.e., there are no impurities inside the wire, and there is no 
scattering from phonons or from the contacts between the wire and the reservoirs at its 
two ends), the zero-temperature conductance is given by Co = 2e^ /h for infinitesimal bias. 
If there are impurities inside the wire which scatter the electrons, then the conductance is 
reduced from Cq. 

A powerful calculational method for studying electronic transport is the non-equilibrium 
Green's function (NEGF) formalism-*^. The advantage of this method is that it treats the 
infinitely extended reservoirs (leads) in an exact way. The derivation of this formalism has 
been based on the Keldysh techniques'^®. Recently a simple derivation of the NEGF results 
based on a direct solution of the equations of motion for a non-interacting system of electrons 
was given in Ref. 10. This method, based on writing quantum Langevin equations, was first 
applied in the case of oscillator systems by Ford, Kac and Mazur^^. It has recently been 
applied in the context of transport^^"^^. 

In the present paper wc point out a particular problem that arises while using the NEGF 
formalism in a situation where there are bound states and there are no additional mechanisms 
for equilibration (such as electron-phonon scattering). We define bound states as states 
whose wave functions decay exponentially as one goes deep into any of the reservoirs. Their 
energy levels lie outside the energy band of all the reservoirs. One expects that the NEGF 
results should reduce to the usual equilibrium results if all the reservoirs are kept at the 
same chemical potential and temperature. This is easily shown to be true in the absence of 
bound states. In the presence of bound states, we show that while the contribution of the 
bound states to the equilibrium density matrix can be obtained within NEGF, the procedure 
is subtle and somewhat ad hoc. It is not clear in this formalism what the mechanism for 
equilibration of the bound states is. Moreover, if bound states are present, the density 
matrix is not unique in the non-equilibrium case. Here we show that the equation of motion 
approach can be used to obtain a clearer understanding of this problem of equilibration 
of bound states and the non- uniqueness of the non-equilibrium steady state. The central 
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results of this paper are as follows. 

(1) We give a simple and general derivation of the NEGF results by the equation of 
motion method for a system without interactions. This is obtained in two different ways: (i) 
from the steady state solution of the equations of motion, and (ii) from the general solution 
involving initial conditions. 

(2) We show that, in the presence of bound states, the exact solution of the wire plus 
reservoir equations of motion (without any additional sources of equilibration) leads to 
steady states which depend on the initial conditions of the wire. 

(3) We show that introducing additional broad-band auxiliary reservoirs (which are very 
weakly coupled to the wire) solves the problem of initial condition dependence. We obtain 
the non-equilibrium steady state properties in the limit where the coupling strength of the 
auxiliary reservoirs goes to zero. We find that the equilibrium density matrix is then unique 
and independent of the properties of the auxiliary reservoirs. But the non-equilibrium 
density matrix depends on the details of the auxiliary reservoirs and on the way in which 
their couplings (to the wire) are taken to zero. 

Bound states have recently been studied in the context of the NEGF formalism^^'^^, but 
to our knowledge, this particular problem of equilibration has not been addressed earlier. In 
this paper, we deal with electronic transport in non-interacting systems modeled by tight- 
binding Hamiltonians. For simplicity, we only consider spinless fermions here although it is 
quite straightforward to include spin. 

The paper is organized as follows. In Sec. II, we discuss the NEGF formalism and present 
the expressions for the density matrix in the wire and the current. In Sec. Ill, we present 
a derivation of the NEGF results using the equation of motion method. This derivation 
is similar to that of Ref. 10 but is a simplified and more generalized version. Starting 
from the full Heisenberg equations of motion of the wire and reservoirs, we derive effective 
quantum Langevin equations for the wire. These equations are solved by Fourier transforms 
to give the steady state solution which leads to expressions for the density matrix and the 
current which are identical to the results obtained from NEGF. In Sec. IV, wc point out the 
problem of equilibration of bound states. In Sec. V, we consider the general solution of the 
equations of motion, as opposed to the steady state solution obtained in Sec. 111. This lets us 
understand better the problem of equilibration in the presence of bound states. The question 
of the approach to the steady state (both in cases with or without bound states) can also be 
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addressed in this approach. In Sec. VI, we describe our method of resolving the problem of 
equilibration of bound states. Namely, we introduce auxiliary reservoirs which are weakly 
coupled to the wire in such a way that the bound states which were earlier localized near 
the wire now extend infinitely into these new reservoirs, and the original bound state energy 
levels now lie within the energy band of the auxiliary reservoirs. The steady state properties 
are obtained in the limit in which the coupling of auxiliary reservoirs to the wire is taken 
to zero. In Sec. VII, we present some numerical results, for a system of a wire with a few 
sites coupled to one-dimensional reservoirs, to illustrate some of the analytical results. In 
Sec. VIII, we briefly consider systems of interacting electrons, and explain why a proper 
treatment of bound states is important for computing the current. In Sec. IX, we make 
some concluding remarks. 

II. NON-EQUILIBRIUM GREEN'S FUNCTION FORMALISM 

In this section, we will briefly discuss the NEGF formalism-'^'^"^'^'-'^'^ for a system which 
consists of a wire connected to reservoirs which are maintained at different chemical poten- 
tials or different temperatures. In the NEGF formahsm, both the wire and the reservoirs are 
modeled by microscopic Hamiltonians. We will use a tight-binding model of non-interacting 
electrons which we will now describe. We use the following notation. For lattice sites any- 
where on the system we will use the indices r, s; for sites on the wire {W) we will use the 
integer indices k, ■ ■ ■; for sites on the left reservoir (L) we use the Greek indices a, v; fi- 
nally, for sites on the right reservoir (R) we use the primed Greek indices a', i/ . We consider 
the following Hamiltonian of the full system: 

n^Y^Hrs clCs (1) 

rs 

= Hw + + Hr + Vl + Vr , 

where Hw ^ ^i^rn , = ^ H^,, c^c^ ^ Ur ^ ^^'u' ^l'^'^' ' 

Im Oiv ct'v' 

Vl = E [ cUa + 4q ] , 

la 

Vil = E [ V3 4ca' + V^i^ clci ] , 
la' 

where cj, denote creation and destruction operators satisfying the usual fermionic anti- 
commutation relations. The parts Tlw, Ti-L and Hr denote the Hamiltonians of the isolated 
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wire, left and right reservoirs respectively, while Vl and Vr describe the coupling of the left 
and right reservoirs to the wire. The main results of NEGF are expressions for the steady 
state current and density matrix in the non-equilibrium steady state. To state these results 
we need a few definitions which we now make. Let G^{uj) be the full single particle Green's 
function of the system but defined between sites on the wire only. (See App. A for definitions 
of the various Green's functions that will be used). Thus if the wire has sites then is 
a, N X N matrix. It can be shown (see App. A) that is given by 

where ^ are self-energy terms which basically model the effect of the infinite reservoirs 
on the isolated wire Hamiltonian. (We will work in units in which Planck's constant h — 1). 
The effective wire Hamiltonian is thus -|- -|- E J which in general will be shown to be 
non-Hermitian. The self energies can be written in terms of the isolated reservoir Green's 
functions gt nijjj) and the coupling matrices V'"'^ (App. A). We get 

E+(a;) = gtiu;) V^^ , E+(a;) = g^{u;) V^^ . (3) 

Finally, let us use the following notation for the imaginary parts of the self energies from 
the two reservoirs, 

r.iu) -^[^I-^t]-V' PL , Vr{u) -^^^R-n\-V^ PR V"' , (4) 

where E~(a;) = Vg~{ijj)V\ and p{uj) — — {1/t:) I m[g^{uj)] is the density matrix of an isolated 
reservoir. 

With these definitions, NEGF gives the following expressions for the density matrix and 
current: 

nim = ( C^m{i)ci{t) ) 

/oo 
-oo 

^ = E ( M K^a cUt)Cm{t) - Vma 4(t)c«(t) ] ) 
ma 

/oo 
doo Tr[ TLG+TnG' ] (/(a;, /x^, T^) - /(u;, /x^, T^)) , (5) 
-oo 

where /(a;,/x, T) = l/[e'^('^~^) -|- 1] denotes the Fermi function. 
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III. DERIVATION OF NEGF BY "EQUATION OF MOTION APPROACH": 
QUANTUM LANGEVIN EQUATIONS AND SOLUTION BY FOURIER TRANS- 
FORMS 



In this section we give a derivation of the NEGF results using the equation of motion 
approach which was developed in Ref. 10. The derivation we present here is a simplified 
and generalized version of Ref. 10. This method basically involves writing the full equations 
of motion of the system of wire and reservoirs. The reservoir degrees are eliminated to give 
effective Langevin equations for the wire alone. Finally the Langevin equations, which are 
linear for the case of a non-interacting system, are solved by Fourier transformations to 
obtain the steady state properties. 

The Heisenberg equations of motion for sites on the wire are: 

Cl^ -i J2 Cm - i ^la Ca - i J2 ^la' ^a' , (6) 

m a a' 

and for sites on the reservoirs are: 

y I 
c^, = -^ E ^a'.' c.' - ^ E K.? Q ■ (7) 

v' I 

We solve the reservoir equations by treating them as linear equations with the term con- 
taining ci giving the inhomogeneous part. Using the Green's functions 
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gtit) = -I e-*^^ e{t) = - du 9l{u)e-^^' , 

ZTT J-oo 
„ 1 /-oo 

gl{t) = -i e-^*^" d{t) = - du: 9^{u;)e-^' , (8) 
we get the following solutions for the reservoir equations of motion Eq. (7) (for t > to): 

Ca{t) = ^ E [9t{t-to)U Cuito) + fdt' E [9t{t-t')U Q(0 , (9) 

Ca'{t) = I E [9Rit-to)]a'.' C,'{to) + fdt' Ebfl(^ " K^/ Ci{t') . (10) 

Plugging these into the equations of motion for the wire Eq. (6), we get: 



Q = E HZ Cm - ivi\ - ^ tdt! E {9l{t-n.. ^^1 c^{t!) 
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-i - ij dt' E [g+(t-t')U, KI Ut') , (11) 

where Jyf = ^ E [5^(^-^o)]a. c,(io) , r^f = ^ E C [^^(^ - ^o)]a'.' c,,(io) ■ 

(12) 

We have broken up the reservoir contributions into noise and dissipative parts; it is clear that 
the wire equations now have the structure of quantum Langevin equations. The properties 
of the noise terms 77^, rf- can be obtained from the condition that at time t = to, the 

two reservoirs are isolated and described by grand canonical ensembles at temperatures and 
chemical potentials given by {Tl, and {Tn, /iji) respectively. Thus we find, for the left 
reservoir, 

( vf^'Mit') ) = E Via [gnt - to)U V^, [gtif - to)U ( cUtoMto) ) . (13) 

Let ipq{a) and denote the eigenvectors and eigenvalues of the left reservoir Hamiltonian 
which satisfy the equation 

E ^^(^) = ■ (14) 

m 

The equilibrium correlations are given by: 

(4(io)c.(io)) = ECH V'^(^) KK^^^L,n) . (15) 

Using this and the expansion {gi^)ycr = —iO{t — to)J2q'ipq{^) '4^q*{<^) e"*'^? in Eq. (13), 
we get 

( V''l{t)vinit') ) = E ^/aM E e^'^^*'*'^ /(Aj,/^L,T^) ] , (16) 

au q 

with similar results for the noise from the right reservoir. Now let us take the limits of 
infinite reservoir sizes and let to ~^ —00. On taking the Fourier transforms 

1 fOO 1 poo 

cicu) = — / dt c{t) e-* , fjicu) = — / dt r){t) e^* , (17) 

ZTT J-00 ZTT J -00 

we get from Eq. (11): 

E [^5im - HZ\ Cm{uj) = E ( [^t{^)]lm + [^U^)]lm ) C^(c^) + (^) + V^^) , 
m m 

(18) 

where E+ = V^g+V^'' , E+ = V^g+V^^ . 
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The terms Ej, are the same self energies that appear in the NEGF formahsm and 
effectively change the Hamiltonian, H^, of the isolated wire, to + E2(a;) + Ej(a;). The 
noise correlations can be obtained from Eq. (16) and give: 

( f,f\u)t^{u') ) = (yV^^V^ f{^,y^L,n) 5(00-00') 

= Tt^iu) fiiu,fiL,TL)6iio-u;') , (19) 
where p^oo) = ^^|^^{q) (?) S{oo - A,^) , 

and similarly for the right reservoir. This is a fluctuation-dissipation relation and shows how 
the noise-noise correlations are related to the imaginary part of the self energy. Finally, we 
get the following steady state solutions of the equation of motion: 

/oo 
doo ci{oo) e-'^' , (20) 
-oo 

with ci{u) = Y: G+iM [ 77^(0;) ] , (21) 

m 

G-M = ^ (22) 

For the reservoir variables, we have from Eqs. (9,10) 

E V,^ c«(^) = fj^iu;) + U^ti^)]im GU^) [ tiuj) + r)^(^) ] , (23) 

a mn 

^ V^^, cAu^) = fi^iu;) + U^Uu^)\im GtnH [ tniu^) + f,^{u^) ] ■ (24) 

a' mn 

1. Steady state current and densities 

Current: From a continuity equation, it is easy to see that the net current in the system 
is given by the following expectation value, 

^ = E( M Ci{t) Ut) - VL Ca(t) ] ) = E2 K^a ( 4W C„(t) ) ]. (25) 

vna ma 

From the steady state solution (Eqs. 20-24), we get 

mix iiiuc ^ "-^ 

doo / doo' e^(--- E ( Gj-^iu;) [ fjf^\oo) + fjf{u^) ] 

X m^') + E [ K(^')U GUuj') ( fij^iu') + v^{u') ) ] ) . (26) 

jk 



Let us consider that part of the above expression which depends only on the temperature 
and chemical potential of the right reservoir. This is given by 

/oo /•oo . 

Imjk 

/oo 
du; Y: GrM[^t{^)U GI,{u;) [r^(a;)],, /(u;, T^) 

Imjk 

/oo 
du Tr[G-{u;) E+(^) G+{u;) r^(w)] f{uj,i^R,TR) . 
-oo 

Taking the imaginary part of this and using this in Eq. (25), we find the contribution of the 
right reservoir to the current to be 

/oo 
duj Tt[VlG+VrG-] f{uj, iiR, Tr) . (27) 
-oo 

[In deriving this result we used the identities = G", = E", = T and the fact 

that (Tr[M])* = Tr[M'''] for any matrix M]. It is clear that on adding the contribution of 
the left reservoir, we will get the net current (left-to- right) to be 

/oo 
duj Tr[TLG+TRG-] {f{uj, i^l, Tl) - /(a;, i^r, Tr)) . (28) 
-oo 

Densities: The density matrix is given by the expression 
nim = ( 4(*) ci{t) ) 

/OO /* OO 

du / du' e^('^--')* E G^u) { [r^f {u) + fjf {u)] [fjf{u') + fjf{u')] ) Gt^{u') 
-oo ^-OO 

/oo 
du [ {G+rLG-)im f{oj, i^L, Tl) + {G+rRG-)im f{0J, I^R, Tr) ] . (29) 
-oo 

The results in Eq. (28,29) are identical to the results from NEGF given in Sec. II. For 
IJ'L — IJ'R and Tl — Tr, we get the expected equilibrium value [see App. (B)] provided that 
there are no bound states. 

IV. PROBLEMS WITH BOUND STATES 

A desirable property of the NEGF results is that they should reduce to the standard 
equilibrium results for the case when both reservoirs are at the same chemical potential and 
temperatures, i.e., ior /il — /J'R — /J' and Tl — Tr — T. For this case the current vanishes. 
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which is the correct result expected from equihbrium. For the density, we get from Eq. (5): 
nim = r du [G+icu) ( ri(a;) + r^(a;) ) G-{u;)]im f{u;,pi,T) , (30) 

The expected equihbrium result can be obtained using the grand canonical ensemble, and 
we get (see App. B): 



eg 



/OO POO f 

du (r^ + Vr) G-\irn f{oj, l^,T) + du- [G+G-J^^ /(a;, T) , (32) 
-OO J-oo TT 

where G^iu) — — rrj — "L, , , ^ , , , . (33) 

Note that here we retain the ie factor in the Green's function. The ie factor in Eq. (31) 
can be dropped since, for uj in the range of interest in Eq. (30), the self energies and 
have finite imaginary parts. (Whenever we introduce e in an equation, it is understood 
that it is an infinitesimal positive quantity which has to be taken to zero at the end of the 
calculation). The first piece in Eq. (32) is identical to the NEGF prediction, while it can be 
shown that the second piece vanishes only if there are no bound states. Thus in the absence 
of bound states we verify that the NEGF results reduce to the correct equilibrium results. 

Let us discuss now the case when there are bound states. Here we refer to bound states of 
the full system of wire and reservoirs. The bound states have energies lying outside the range 
of the reservoir levels, and the wave functions corresponding to them decay exponentially as 
one goes deep into the reservoirs. Tlicy are obtained as real solutions of the equation 

[//^ + E+(A6) + E+(A6)] *6 = A6*6 . (34) 

From the form of the Green's function, it is clear that the second term in Eq. (32) is non- 
vanishing whenever the above equation has real solutions, and it can be shown to reduce 
to the form 



n 



eq—bnd 
Im 



= Y: Mmim) f{Xb:Pi:T) . (35) 

6 

Thus there seems to be a problem of equilibration in the NEGF formalism whenever bound 
states are present. This problem can be fixed in the following way. In the NEGF results 
let us add an extra infinitesimal part ea/7r to the matrix ra{ui) corresponding to the a^^ 
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reservoir (a — L, R). Thus the NEGF result for the density matrix is modified to 

/oo 
-oo 



-oo 
roo 



+ r du[^-^ {G+G-)im f{u;,i^L,TL) + - {G+G-)im f{0U,I^R,TR) ], 

J-oo TT TT 

with G~^(u}) — —77-. — ^1 , , ^ , , ■ . (36) 

The second integral gets contributions from bound states. Using the identity 

lim , = lim tt — — — Siuj - A^) , (37) 

we finally get the following contribution from the bound states, 



n 



6 ei,ii-^0 cl + eR 

For f^L = f^R = and = Tr = T, we now get the expected equilibrium result of 
Eq. (35). For the non- equilibrium case, however, the situation is somewhat unsatisfactory 
since the density matrix depends on the ratio cl/^r, and so the bound state contribution 
is ambiguous. More importantly, in this approach it is not at all dear as to what the exact 
physical mechanism for equilibration of the bound states is. 

In the following sections, we will examine this particular question of equilibration of 
bound states more carefully. 



V. GENERAL SOLUTION OF EQUATIONS OF MOTION 

In this section, we will consider the general solution of the equations of motion. Unlike 
the Fourier transform solution obtained in Sec. Ill, the general solution also involves the 
initial conditions of the wire. The advantage of the equations of motion method is that 
it can address issues such as that of approach to the steady state, it can be numerically 
implemented (after truncating the reservoirs to a finite number of sites) , and it does not rely 
on any ad hoc ie prescriptions for the reservoirs. We will show that the expression for the 
density matrix can be derived without any difficulties if there are no bound states. Further, 
the problems which arise if there is a bound state become quite clear in this approach. 

Let us again consider a system in which a wire is connected to two reservoirs, with 

n ^ ^ Hrs 4 c, . (39) 

r,s 
11 



As before we will use the label r, s to denote any site in the system, and j, a and a' to 
denote sites in the wire and left and right reservoirs respectively. Instead of eliminating the 
reservoir degrees of motion and writing Langevin equations for the wire, we will now deal 
with the full Heisenberg equations of motion for the system which is given by 

Cr = - i ^ Hrs Cs . (40) 
s 

At the initial time i = 0, we assume that the two reservoirs are in thermal equilibrium 
at temperatures and chemical potentials T^ r and /iL,R, and the couplings V^'^ are zero. 

As before, let and 7/;^, with a = L, R. denote the eigenvahies and eigenfunctions of the 
isolated reservoir Hamiltonians. At t = 0, the density matrix for sites on the left and right 
reservoirs are given respectively by 

(4(0)c.(0)) = E /(AJ,/.^,T^) , 

(i(o)c.,(o)) = y: V'f /(Af,/x«,r«) . (41) 

q 

We also assume that &t t — 0, there are no correlations between the two reservoirs and 
between the wire and the reservoirs. Thus 

( 4(0) c«(0) ) = ( 4(0) cAO) ) = ( 4(0) c^,(0) ) = . (42) 

Finally, for sites in the wire, we cannot unambiguously assign a value to the density matrix 
(c|(0)c^(0)) at t — 0, since the wire is isolated from everything else at that time. Under 
certain conditions, we will see that in the steady state (defined as the limit t — > 00), the 
density matrix of the wire will turn out to be independent of its value at t = 0. 

Let us now suddenly switch on the couplings Va of the wire to the reservoirs at t = 0. 
For t > 0, the solution of the equations of motion is given in matrix notation by 

c{t) = I g+{t) c(0) , (43) 

where Q^{t) = —ie~'^^*9{t), and c(t) = [ci, C2..., ca^J^, where Ns is the total number of sites 
in entire system. For a point / in the wire, we then have 

ci{t) = E sLit) uo) + E c.(o) + E G^.' c.'(o) . (44) 

m V n' 
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Hence the density matrix in the wire is given by 

nimit) = {cl{t)ci{t)) 

= E o^jit) (c|(o)c,(o)) gr^it) + E g^^it) (ct(o)c«(o)) 

ij OLV 

+ E ^;a'(^) (4'(0)Ca'(0)) ^.-^(t) . (45) 
a'v' 

Now let ^E'Q(r) and Aq denote the eigenfunctions and eigenvalues respectively of the system 
of coupled wire and reservoirs. Thus 

Hrs ^q{s) = AQ*Q(r) . (46) 

s 

Then the full Green's function can be written in the following way (for t > 0): 

Q 

= -z ^ ^b{T)n{s) e-'^^' - t [du; p'i.M e— * , (47) 

b ■' 

where Af, refer to bound states, and the density matrix p'^ is given by a sum over the 
extended (continuum) states of the system ^^^(ci;) = X)q ^q(^)^q(s)<5(w — Aq). In the limit 
i — > oo, the second term vanishes (this follows from the Riemann-Lebesgue lemma, see 
Ref. 19), and we get a contribution only from bound states. Thus 

hm gUt) = -^ E "^birWbis) e"^^^* . (48) 

t— >0O 



Thus if there are no bound states in the fully coupled system then, for any two sites of the 
system, ^^(t) vanishes as t oo. From this it follows that the contribution of any finite 
number of terms in Eq. (45) vanishes in the steady state. Since the wire consists of a finite 
number of sites, this means that the initial density matrix of the wire (ct(O)cj(O)) will have 
no effect on the steady state density matrix nim{t — > oo). The reason that Eq. (45) does 
not vanish in the steady state is that it gets contributions from the reservoirs which have an 
infinite number of sites. On the other hand the situation is quite different if there arc one or 
more bound states in the problem. In that case, clearly, individual matrix elements of g^ 
may not vanish in the steady state, and the initial state of the wire makes a contribution to 
the long time density matrix. 

We now look at the contribution of the reservoirs in Eq. (45). Dropping the subscript a 
for the moment, let us look at the contribution of any one of the reservoirs. Using Eq. (41), 
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this can be written as 



^res 



/oo 
-oo 



Im 

av 

oo 



did f{uj,ij,,T) Iimiuj,t) , (49) 

'J 

where 7,^(u;,i) = Y.^^ i E i^ii^m^^) S{u; - X,) ] g-^{t) . 

av q 

We now transform to the frequency dependent Green's function using the relation 
Q'^it) — /f^ |^^+(a;)e~"^*. As shown in App. A, we can express the Green's function 
elements Qji,{oj) in terms of Gfj{uj) and the isolated reservoir Green's functions g'^^{uj). We 
get 



and g-i{uj) = Y 9aM 4 Gjiiuj) . (50) 
We will also use the following result which follows from an eigenfunction expansion of g~^{u!), 

a a q' ^ + ^^- \' UJ + 16 - Ag 

Using Eqs. (50-51), we can write Iim{ui,t) in Eq. (49) as 



CO" — le — u 



27r 7 ^ a;' + ie - a; ^ ^ 27r y a;" - ie - a; ^ ^ 

The first integral in the above equation can be evaluated as follows: 

J uj' -\-ie — uj ^ ^ 27r 



2t: J cu' + ie-u ^ ' 2tt J-oo {uj' + ie - uj){uj' + ie' - Xq) 



Q 



= -^GJ(a;)e-* + *q(0^q(j) ^ ^ _ ■ (53) 



In the limit t — > cxd, only the bound states contribute to the summation (over Q) in the 
expression above. Hence we get 



lim 

t— »oo 



ZTT J CO + le — CO ■' ^ U — Ah 
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Note that we have dropped the ie factor in the denominator of the second term. This is 
because of the following reason. In the expression Eq. (52), the presence of the density of 
states of the reservoirs p{uj) means that we will be interested only in values of cu lying in 
the continuum of reservoir eigenvalues. Using Eq. (54), we finally get the following result 
for the contribution from the a*^ reservoir to the density matrix in the long time limit: 



n 



Im 



/oo 
duj [ G+{u;)T-{u;)G-{u;) ^ /(c^,/ia,T„) 
-oo 

bb' ji {uj - Xb)K^ - M') 

(55) 



where 6, h' denote bound states, and we have again dropped terms in which the integrands 
contain time-dependent oscillatory factors. 

In the absence of bound states, we see that the contribution of the reservoirs is identical 
to that given by NEGF and also by the Fourier-transform solution in Sec. III. However, 
in the presence of bound states, there is an extra contribution which is not obtained in 
the other methods. Importantly, this extra part is, in general, time-dependent so that the 
system never reaches a stationary state. Also, we do not recover the equilibrium results for 
the case when all reservoirs are kept at equal temperatures and chemical potentials. The 
reason for this is that the reservoirs are not able to equilibrate bound state energy levels 
since these levels lie outside the range of the reservoir band of energies. 



VI. A SOLUTION TO THE BOUND STATE PROBLEM 

One way to solve the bound state equilibration problem is to introduce two auxiliary 
reservoirs L and R\ these must have the same temperature and chemical potentials as the 
original reservoirs L and i?, but they must have a large enough bandwidth so that any 
bound states of the original Hamiltonian H. lie within that bandwidth. (Briefly, the idea is 
that since the bound states of H. are no longer bound states in the expanded system, the 
arguments given in Sec. IV imply that the density matrix of the wire will no longer depend 
on the initial density matrix). The two auxiliary reservoirs must be coupled very weakly to 
the wire, so that they do not greatly alter the energies of Ti and the corresponding wave 
functions within the wire and the original reservoirs. We will eventually take the limit of 
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that coupling going to zero. The main purpose of the auxihary reservoirs is to equihbrate 
any bound states which H might have. Let us now see how all this works. 
We consider a new Hamiltonian of the form 

^new + Hi + Hr + Vi + Vn, (56) 

where H ^ Hw + nL + nR + VL + VR . 

Ha denotes the Hamiltonian of the auxiliary reservoirs (a = L,R), and denotes the 
coupling of those reservoirs to the wire. The auxiliary reservoirs will also be taken to be 
lattice systems with tight-binding Hamiltonians, but with a hopping amplitude which is 
sufficiently large. We will take the auxiliary reservoir L to have the same temperature and 
chemical potential Tl and fii as the reservoir L, and similarly for the reservoirs R and R. 

Following the derivation of Sec. Ill, we find that the full Green's function on the wire is 
now given by 

^^^^^ ^ u; - - Etiuj) - E+(^) - E±(c.) - Et(^) ' ^^^^ 

With the auxiliary reservoirs present, there are no longer any bound states, and the density 
matrix in the wire can be written as 



/oo _ _ 

d(jj [G+(rL + Ti)G~]im f{uj, hl, Tl) 
-oo 

/oo _ _ 

du [G+{VR + VR)G-\imf{uj,iJiR,TR) 
-oo 



We want to eventually take the limit in which the couplings of the wire to the auxiliary 
reservoirs Va go to zero. We will therefore treat these couplings perturbatively. Let us 
first break up the above integral over uj into two parts, one with uj going over the range of 
the original reservoir band, and the other containing the remaining part over the range of 
the auxiliary reservoir. We assume that the original reservoir bandwidths are in the range 
[£^1,-^2] while that of the auxiliary reservoirs are [Ei,E'2\; we choose the latter such that it 
contains the original bandwidth, i.e., Ei < Ei and E2 > E2. Thus we write: 



n]^ = I duj[G+{TL + Ti)G-]imf{uj,iiL,TL)+ I ' duj [G+{T r + T R)G-]imf{uj, iir,Tr), 

JEi JEi 

nL = rdu;[G+riG-]imfiuj,fiL,TL) + ^ duj [G+V RG'^f {uj , ^.Tr) 

JEi JEi 

+ r du[G'-T-LG-]imf{u,iiL,TL) + r du[G+VRG-]imf{u,iiR,TR) . (58) 

J E2 J E2 
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Now, in the first part {-n}) we note that both and have finite imaginary parts. Hence 
on taking the hmit V^.J? ~^ 0; we immediately get Tij^ — > and = Hence 

nL = rdij[G+rLG-]imf{iu,fiL,TL)+ rduj[G+riiG-]imf{iJ,fiR,Tit), (59) 

JEi JEi 

which is just the usual reservoir contribution. 

We will now show that the second part contributes to bound states only. As before, 
let us denote the eigenvalues and eigenf unctions of the combined system of the wire and 
reservoirs (but not the auxiliary reservoirs) by Xq and ^q, where the label Q can have both 
a continuous and a discrete part (corresponding to bound states) . We will assume that the 
auxiliary reservoirs have been chosen such that all the values of Xq lie within the bandwidth 
of the auxiliary reservoirs. The coupling to the auxiliary reservoirs cause the energy levels 
to shift from Xq to Xq + 6Xq. To first order in Vg,, we find that 

SXq = (*Q|Sj(AQ)+Et(AQ)|*Q) = Y: nU) Pj + Et],,*Q(fc) . (60) 

jk 

The imaginary part of SXq is therefore given by 

Im SXq = -tt^ [Ti + T^],,, ^q(A;) . (61) 



Thus for the auxiliary reservoir a we get 

km km Q ^-Xq-SXq q, OJ-Xq'-SXq, 



(62) 



Since wc arc considering values of u outside the range [Ei, E2], it is clear that, in the limit 
Fa —> 0, only bound states with Xq — Xqi — will contribute. The above then gives 

■-■-» t - V (- - A* - Re S\f + (Im SX,r 



Wrl + r^k) (63) 

where we have used the fact that both real and imaginary parts of ^A^ are small, with the 
imaginary part given by Eq. (61). 

Putting everything together, we finally find that the total density matrix in the wire due 

to all four reservoirs is given by 

nim = J^' dw [ ( G+{u;)rL{u;)G-{u;) )im /l(c^) + ( G+ {u;)r r{u;)G- {co) )im fniuj) ] 
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where the sum over b runs over all the bound states of H. In the equilibrium case [i.e., for 
fL{^b) — fRi^b)] Eq. (64) shows that the contribution of a bound state b to the density 
matrix is independent of details of the auxiliary reservoirs such as the quantities (^5|rj^|\l/ft) 
and (^6|r^|\l't). But in the non- equilibrium case [i.e., for fii^b) 7^ fni^b)], the contribution 
of a bound state does depend on details of the auxiliary reservoirs, namely, on the ratio 
(^b|r£|^;,)/(^(,|r^|^;,). This is perhaps not surprising. Some quantities in equilibrium 
statistical mechanics may be independent of the mechanism for equilibration, while the 
same quantities in non-equilibrium statistical mechanics may depend on the details of that 
mechanism. 

The similarity between Eq. (38) and the bound state contribution in Eq. (64) implies that 
we can think of the auxiliary reservoirs as providing a justification for the prescription 
which was introduced in the NEGF formalism in Sec. II. 

We note that bound states do not carry current, and so the expressions for current remain 
unchanged. Later we will point out that in the presence of electron-electron interactions, 
even the current is likely to be affected by bound states. 

VII. NUMERICAL RESULTS 

It is instructive to consider the following simple examples where we can explicitly see the 
problem of equilibration of bound states and its resolution by the introduction of auxiliary 
reservoirs. 

A. Wire with a single site 

We will first consider a system where the wire (W) consists of a single site with an on-site 
potential and connected to one reservoir (i?) and one auxiliary reservoir [R). The reservoirs 
are one-dimensional semi-infinite electronic lattices. The full Hamiltonian is thus given by 

H = Hw + 'Hr + Hr + Vr + Vr , 
Hw = VcIcq , 

-2 

a=—oo 

Vr = -7' (cLiCo + 4c_i) , 
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a=—oo 

Vr = -f(cljCo + 4c_i) . (65) 



Let us first give the exact equilibrium results for the system of wire and the single reservoir. 
The reservoir will be assumed to have a chemical potential /i and temperature T. Prom 
App. B, we get for the density at the single site on the wire 

where = f^du 7" |G+o(cu)P p{uj) /(cu,/x,T) , 

/'OO c 
ciu;-|G+oH|V(^,/^,r), (66) 
-00 TT 

and GnJuj) = — — ——- . (67) 

We need g'^i^^) and piuj) which are the reservoir Greens function and the reservoir density of 
states respectively, both evaluated at the site a = —\. The eigenvalues and eigenfunctions 
of the reservoir Hamiltonian are given by 

\q = —27 cosg , and i^qipi) = \/2sin(g a) , (68) 

where q lies in the range [0, vr], and n = —1, —2, ■ ■ ■. The wave functions are normalized so 
that {'4'q\'4'q') = T^5{q — q'). Hence we get for the required reservoir Green's function, 

gV) = for M<27 

= for |a;| > 27 . (69) 
Similarly, for the auxiliary reservoir we get. 



-9^{^) = T.[^-^{^-^f''] for H < 27 



7^27 ^ 472 

i[^ + (^-l)V2] for c.<-27 
7 27 472 
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= for |a;| > 27 . (70) 

Bound states: We choose the parameter values 7 = 7' = 1. For the coupled system 
of wire and reservoir, we again get a continuum of states identical to the original reservoir 
levels. In addition, for y < — 1, we get a bound state whose energy is given by = V + l/V, 
and the wave function at the O*'' site is ^^(O) = (1 — 

System without auxiliary reservoir: In our numerical example we take V — —2, fi — 
0, and P — 2. Thus, from Eq. (67), the expected equihbrium results are 

<' = f ^ ,"'"1/'"-':^' = 0.1641... , 

nf = *;.(0)V(A(„//,r) = 0.7449... , 

rfi = n^^ + Ufe" = 0.9091... . (71) 

We have numerically solved the equations of motion for a reservoir with L = 500 sites. 
In Fig. 1, we show the contribution of the reservoir to the density at the 0*^ site as a 
function of time by a solid line. We find that it quickly settles at a value of about ~ 
0.314. The contribution coming from an initial density equal to one at site 0, namely, 
lim^^oo GQo(t)GQo(^) Eq. (45), is shown by a dashed line. This does not vanish with time 
but settles at a value of about n}"^ — 0.563. This means that the steady state density on 
the wire depends on the initial density. Using the results of Sec. V we can understand the 
different contributions to the density. From Eq. (48), we get the contribution from the initial 
density as 

n^" = |*b(0)|^ = 0.5625 , (72) 

which is close to the numerically obtained result. From Eq. (55), we see that there arc two 
parts to the contribution from the reservoir levels. The first part is a contribution to the 
density at site arising from the extended states; this is identical to the equilibrium result 
— if^ — 0.1641.... The second part is from the bound state; from Eq. (55) this is given 

by 



nt. = 1*6(0)1^ r du;-^^^f{u;,pi,T) = 0.1500... 

J-2 [00 — \h)^ 



'■bnd 



{u; - XbY 
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Hence we get n = nf^^ + = 0.3141.... We note in Fig. 1 that for large times which are 
of the order of the reservoir size, both the contributions start deviating from their steady 
values; we will comment more on this below. 
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FIG. 1: Time evolution of the contributions of a single reservoir (n^) and the initial condition on 
the wire (n^") to the density at that site, for 7 = 7' = 1, /3 = 2, /i = 0, and V = —2. 
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FIG. 2: Time evolution of the contributions of the two reservoirs (n^) and the initial conditions on 
wire (n^") to the density at that site, for 7 = 7' = 1, 7 = 2, 7' = 0.2, /? = 2, /i = 0, and V = —2. 

System with auxiliary reservoir: We now introduce an auxiliary reservoir with the 
same temperature, chemical potential and number of sites as the original reservoir. However, 
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we take the hopping amphtude in this reservoir to be 7 = 2, so that its bandwidth of [—4, 4] 
includes the energy of the bound state of the original system. We take the coupling of the 
auxiliary reservoir to site to be 7' = 0.2. We now solve the equations of motion of this new 
system consisting of two reservoirs and one site. Fig. 2 shows the total contribution of the 
two reservoirs by a solid line, and the contribution coming from an initial density of one at 
site 1 by a dashed line. We see that the solid line approaches a value of 0.907 (this deviates 
slightly from the equihbrium value of 0.909 because 7', though small, is not zero), while the 
dashed line vanishes around the same time. Using the results in Sec. VI, we can compute 
the contributions of the reservoir and the auxiliary reservoir to the net density. We get 



Hence the total density is n — n + n — 0.9072... which is consistent with the value 
obtained from the numerics. Comparing with Eq. (71) we note that the auxiliary reservoir 
only contributes to equilibration of the bound state. 

In Fig. 3, we show an example where the potential V = —0.5 is such that there are no 
bound states. In this case we find, as expected, that the wire equilibrates even with a single 
reservoir. Adding an auxiliary reservoir leaves the density essentially unchanged. 

The time scale of approaching these steady values is of the order of r = 7/7'^ = 50. This 
is due to the fact that the self-energy of the auxiliary reservoir is proportional to 7'^/7. This 
governs the rate at which the auxiliary reservoir fills up site 1, and it is inversely proportional 
to r. For times t >> r, we obtain the correct equilibrium value of the density at site (up 
to a small error due to the finiteness of 7'). Further, the vanishing of the dashed line means 
that the steady state density at site does not depend on the initial density at that site. 
Once again, we see in Fig. 2 that for large times of the order of the reservoir sizes, both the 
contributions start deviating from their steady values. 

The deviations which begin appearing in Figs. 1-3 at large times of the order of the 
reservoir sizes can be understood as follows. After the wire is suddenly connected to one 
end of a reservoir at time t — 0, there is a recurrence time of the reservoir which is given 
by the time it takes for the effect of that connection to propagate to the other end of the 
reservoir, and then return to the wire. If the reservoir has a size L, and the Fermi velocity 
in the reservoir is vp (this is equal to 27sinA;F and 27 sin /c^^ for the original and auxiliary 




(73) 
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FIG. 3: Time evolution of the contributions of a single reservoir (n^) and the initial conditions 
on wire {n^"') to the density at that site, for 7 = 7' = 1, /3 = 2, /x = 0, and V = —0.5. The inset 
shows the time evolution for the case with two reservoirs. 

reservoirs respectively), the recurrence time is given by 2L/vf- Since we have taken fi = 
(half-filling), we have kp = it/2, and vp = L/'j and L/7 in the two reservoirs respectively. 
In Figs. 1-3, we see deviations occurring at a time given by L = 500 due to the original 
reservoir, while in Fig. 2, we also see deviations occurring at a time equal to L/7 = 250 due 
to the auxiliary reservoir. 

The presence of the recurrence times imply that one must take the limits of certain 
quantities going to infinity in a particular order, in order to numerically obtain the correct 
steady state values of the density matrix. The coupling 7' of the auxiliary reservoirs to the 
wire must be taken to zero so that the steady state values do not depend on the details of 
those reservoirs (at least in the equilibrium case); this means that r must go to infinity. On 
the other hand, the recurrence time 2L/vf must also go to infinity. The steady state will 
then exist for times t which satisfy r << t << 2L/vp. 

Approach to equilibrium: Let us briefly discuss the way in which equilibrium is 
approached at times which are small compared to the recurrence times discussed above. 
One again, there are two different time scales here, one for the original reservoirs and the 
other for the auxiliary reservoirs (if there is a bound state present). For the original reser- 
voirs, let us consider an integral of the form given in the second term in Eq. (47), namely, 
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duj p(oi))e~'^^, where the function has no singularities in the range [Ei,E2\. Let 
us also suppose that p{u;) has the power-law forms (a; — -Bi)"^ and (£'2 — a;)"^ at the two 
ends of the integral, where cti, 02 > —1 so that the integral exists. Then one can use 
the method of steepest descent to show that in the limit t 00, the above integral only 
gets a contribution from the end points, and those contributions vanish as e-^-^iy^^i+i and 
^-iE2t i^a2+i respectively^^. For the original reservoirs with 7 = 7' = !, this approach to 
equilibrium occurs at times of the order I/7, and it is therefore hard to see the power-law 
fall-off in Figs. 1 and 4. For the auxihary reservoirs, the time scale of equihbration is given 
by a different expression if there is a bound state present. The relevant integral is then 
given by /J|^ du p{uj)e^'^'^^, where the function p{uj) is large in the vicinity of the bound state 
energy A5, namely, p{uj) ~ {^/"^ / / [{uj — A^)^ + (7'^/7)^] where 7/7'^ is large. We then see 

that the dominant contribution to the integral comes from the vicinity of a; = and is 

-'2 /- 

given by e~* ''' We can see this exponential decay in Fig. 2, till the effect of the shortest 
recurrence times starts becoming visible. 

B. Wire with two sites 

In the presence of two bound states, the steady state properties can be time dependent. 
We now illustrate this with the example of a two-site wire. We take both sites (labeled 
and 1) to have an on-site potential V — —3. The site is connected to a one-dimensional 
semi-infinite reservoir (going from —1 to —L as in the previous example), while site 1 is only 
coupled to site 0. The Green's function is thus 



Let us look at the density at site in the long time limit. From the equation of motion 
solution in Sec. V, we expect the density to have contributions from both the initial density 
matrix on the wire (n^") and from the reservoir (n^). We choose the initial wire density 
matrix to be diagonal, with (cqCo) = 1. Then from Eq. (47) we get 




(74) 



n 



In 



*6i(0)|" + 1*6,(0)1" + 2 |*6i(0)n*6,(0)|2 cos{X,,-X,,)t . 



(75) 
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The contribution from the reservoirs consists of two parts: one corresponding to the extended 
states n^j, and the other to the bound states n^^. These are given by 

r2 



I, = J_^dcj \G+,\' Pico) f{u;,i,,T) 

R _ i.T, inM'i J, . Pi^)fi^^(^^T) 4 p{u;)f{u;,n,T) 



J-2 [UJ — Abij J-2 [UJ — Ab2) 

+ 2 \%M\'\^>^M\' cos (A,, - A,Jt r du; ^^""Y^!^/ , ■ (76) 

J-2 [UJ - Xb^){uj - Xb^) 
The bound state eigenvalues are obtained by solving, for u; — —2 cosh a < 0, the equation 

(tu ~ V — q^(uj) 1 \ 
] = 0. (77) 
1 uj-v J 

With g'^{u;) — — this gives the equation 

e^" + 2ve'^" + -y^e" + v = , 

which, for V = —3, has two solutions with a > 0. These give the eigenvalues A^^ = —2.257... 

and A;,2 = —4.137.... The eigenvectors can also be found easily after normalizing them over 

the entire system. At site 0, we get 

^ (l-e-^")(t; + 2cosha)^ 
' ^' (l-e-2«) + (^; + 2cosha)2 " 

This gives I^^J^ ^ 0.2947... and l^'b^P = 0.5421.... Then from Eq. (75) and Eq. (76) we get 

n^" = 0.381 + 0.3196 cos(1.879t) , 
nf,, = 0.1097 , 

= 0.05654 + 0.045 cos(1.879i) . (78) 

We verify these results by an exact numerical solution of the time evolution of the system 
with the two-site wire and a one-dimensional reservoir with L — 600 sites. In Fig. 4, we 
compare the analytic results with the numerical solution and find very good agreement 
between the two. The inset shows the long time behavior; we sec that the effect of the finite 
size of the reservoir shows up at the recurrence time 2L/vf = 600. In Fig. 5, we show the 
effect of adding a weakly coupled (7' — 0.3) auxiliary reservoir with 7 = 3. As expected, 
in this case there is no contribution from the initial density of the wire, and the reservoir 
contribution gives the equilibrium value (till the effect of the recurrence time of the auxiliary 
reservoir, L/^ — 200, shows up). 
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FIG. 4: Time evolution of the contributions of the reservoir (n ) and initial conditions on wire 
(n-^") to the density at the first site of a two-site system, for 7 = 7' = 1, /3 = 2, // = and v = —3. 
The analytic predictions are also plotted. 




FIG. 5: Time evolution of the contributions of the two reservoirs (n ) and the initial density on 
site (n^") to the density at the first site of a two-site system, for 7 = 7' = 1, 7 = 3, 7' = 0.3, /3 = 2, 
// = 0, and V = —3. 

VIII. INTERACTING SYSTEMS 

Let us briefly consider how interactions can be studied within the NEGF formalism. 
We note that this formalism works with a one-particle Hamiltonian, e.g., the self-energy in 
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Eq. (2) is given in terms of the energy of a single electron which is entering or leaving the wire. 
One way to deal with interactions is therefore to do a Hartree-Fock (HF) decomposition. 
For instance, if we have a Hubbard model with an on-site interaction between spin-up and 
spin-down electrons of the form Ucl^^Cn^cl^^Cni, we can approximate it by as 

U [ (cltCnT) ci^Cni + cl^Cn^ (cl;C„;) - (4^C„t) (c1|C„|) ] . (79) 

We see that the first two terms modify the on-site potential. A self-consistent NEGF calcu- 
lation can then be implemented as follows^°. 

• Start with the Hamiltonian with no interactions, and calculate the density matrix. 
The diagonal elements of the density matrix give the densities at different sites. 

• Use the HF approximation to compute the Hamiltonian with interactions, and use 
that to calculate the density matrix again. 

• Repeat the previous step till the density matrix stops changing. 

• Use the converged density matrix to calculate the site densities or the current. 

The important point is that in an interacting system, the density affects the on-site 
potential and therefore the current. If there is a bound state present, it will affect the 
current by modifying the density, even though the bound state does not directly contribute 
to the current. Hence a proper treatment of bound states is necessary in order to calculate 
the current. Numerically, it is found that the presence of a bound state has a significant 
effect on the conductance of an interacting system^'^. 

The effect of on-site Coulomb interactions on the conductance through a quantum dot has 
been studied earlier in Ref. 21 using a self-consistent truncation of the equations of motion 
for the Green's function. The density on the dot plays an important role in that analysis; it 
is therefore clear that the presence of a bound state would affect the conductance. 

IX. DISCUSSION 

In conclusion, we have pointed out that in the presence of bound states and no additional 
mechanisms for equilibration, the NEGF formalism gives a non-unique density matrix. We 
have shown that the equation of motion approach, which gives a simple and straightforward 
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method of deriving the results of NEGF for non-interacting systems, also provides a clear 
understanding of the ambiguity in the density matrix when bound states are present. 

We have then presented a way of resolving the ambiguity which arises when there are 
bound states. Namely, for each reservoir, we introduce an auxiliary reservoir with the 
same temperature and chemical potential; the auxiliary reservoir is taken to have a larger 
bandwidth (so that the bound state energy lies within that bandwidth) , and a much smaller 
coupling to the wire than the original reservoir. We then find that the contribution of the 
bound state is completely determined by the properties of the auxiliary reservoirs. If we let 
the couplings of the auxiliary reservoirs tend to zero, we get the correct equilibrium density 
matrix including the contribution from the bound states. In particular, we find that the 
bound state contribution is unique and independent of the details of the auxiliary reservoirs 
in the equilibrium case. However we again find that in the non-equilibrium case, the bound 
state contribution is not unique and depends on the details of the auxiliary reservoirs and 
the way in which the limit of zero coupling is taken. 

One might suspect that the sudden switching on of the wire-reservoir interactions could 
be a reason for the problem of equilibration. However we have also studied the case where 
the coupling is switched on adiabatically and verified that the equilibration problem remains. 

For non-interacting systems, bound states do not directly contribute to the current, and 
so transport properties are not affected. However in the presence of electron-electron inter- 
actions, a simple mean-field treatment suggests that bound states affect the local density 
which in turn affects the current. Thus transport properties of interacting systems can be 
affected in a non-trivial way in the presence of bound states. 
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APPENDIX A: GREEN'S FUNCTION PROPERTIES 



The single particle Green's function for the coupled system of wire and reservoirs is given 



by 



It satisfies the equation of motion 



.dg 



+ 



dt 



Hg+ = s{t) . 



The Fourier transform g~^{u;) — /f^ dt ^"'"(i)e*'^* is thus given by 

1 



which can also be represented as 



{cu + ie)i - H ' 



(Al) 



(A2) 



(A3) 



Q 



uj + ie- Xq 
OJ - \o 



ITT 



where and Xq respectively denote eigenfunctions and eigenvalues of the full system. Now 
what we need is the part of the full system Green's function defined between points on the 
wire. Let us try to express this part, which we will genote by in terms of the isolated 

reservoir Green's functions given by: 

1 , . . 1 



(A4) 



of the isolated reservoirs. We rewrite the Green's function equation, Eq. (A3), by breaking 
it up into various parts corresponding to the wire and reservoirs. Thus we get 

({uj + ie)i-H^ -V^ \ ( G+ G+r GtrA //0 0^ 



V 



Lt 



\GrW Gj^j^ Gj^ J 



{u + ie)I -H^ 
- V^^ (00 + ie)i -H^ J 

This gives the following equations: 

[{uj + te)i-H'^]G+ - V'^Gtw - V^G^^ = /, 

- G+{uj) + [ (u; + ie)i - H"^ ] G+w 

- V""' G+iuj) + [{uj + ie)i - //^ ] 



0/0 

Vo /; 



0, 
. 



(A5) 
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Solving the last two equations gives G^|^(a;) = gti^) ^'^i^) ^mvi^) 

9r{^) G'^{u;). Using this in the first equation then gives 



1 

u + ie- -El 
where E+(a;) = V^g+iuj) V^' , E+(c.) = gl{uj) V"' 



GH^) - rsw ^+ ' (A6) 



APPENDIX B: EQUILIBRIUM PROPERTIES 

Let us calculate the expectation value of the density matrix for the case where the entire 
system of wires and reservoirs are described by a grand canonical ensemble at chemical 
potential /i and temperature T. For points on the wire this is given by 

nZ = ( 4 Q )e. = YMQWrniQ) f{^Q,f^,T) (Bl) 

Q 

/oo 
Q 

/oo 1 
du; —[G-{u;) - G+{u;)U f{uJ,ii,T) . (B2) 
-oo 2m 

Now from Eq. (A6) we get 

- (G-)-i = (E;^ - E+) + (E], - E+) + 2te = 2'k%{Y + r«) + 2ie , 
^ G~ -G+ ^ 27ri G+{rL + Vr) G' + 2ie G^G' . (B3) 

Therefore we get for the equilibrium density matrix 

/oo fOG f 

du [G+iFL + rR)G-]im ficu, i^,T)+ du - [G+G-]im fiou, T) . (B4) 
-co J —CO TT 

The second part is non- vanishing only if the equation 

Det [ui -H^ - i:l{uj) - E^(a;)] = (B5) 

has solutions for real uj. These correspond to the hound states of the coupled system. As 
usual, the limit e ^ is implied in Eq. (B4). The second term in Eq. (B4) survives in that 
hmit if there is a bound state, since e G^G~ ~ e/[(a; — A^)^ + e^] — t^5{ijj — A^) in that case. 
In fact it is easy to see that the second term in Eq. (B4) reduces precisely to the following 
form: 

dw - [G+G-]i^ f{uJ,fx,T) = Y "^libWmib) f{K,f^,T) , (B6) 

-oo TT ^ 
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which can also be seen directly from Eq. (Bl). 
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